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摘 E: 采用 Morlet 小 波 分 析 法 对 2019 年 乌 梁 素 海 冰 下 水 体 的 温度 ,溶解 氧 (DO ) 浓 度 进 行 趋势 分 析 。 结 果 表 明 
CL) 湖泊 水 体温 度 存在 分 层 现象 ,同时 ,DO 浓度 也 出 现 不 同 程度 的 分 层 现象 。 表 层 DO 浓 度 近 一 个 月 时 间 内 保持 不 
变 ,中 层 、 底 层 DO 浓 度 随 冰 封 时 间 的 增长 开始 出 现下 降 , 中 层 下 降 速 率 为 0.32 mg L d ,底层 下 降 速 率 为 0.05 mg 
Ltd! (2) 表层 水 温 主 周期 是 25 d、16 d, 在 第 一 主 周期 下 ,周期 约 为 15 d;DO 浓 度 主 周期 是 23 d, 平 均 周 期 为 15 d; 
中 层 水 温 主 周期 是 14 d ,DO 浓度 主 周期 是 19 d、7.5 d, 在 第 一 主 周期 下 ,水 温 平 均 周期 约 为 9d,DO 浓 度 平均 周期 为 
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5 d; 底 层 水 温 主 周期 是 14 d、2.5 d, DO 浓度 主 周期 是 24 d 10 d。 水 温 变 化 周期 约 为 8.5 d, DO 浓度 


平均 周期 7d 左 


右 。(3) 小 波 相干 谱 结 果 表 明 :湖泊 水 体温 度 与 DO 浓度 存在 复杂 关系 ,表层 、 底 层 水 温 与 DO 浓度 呈 反 相位 关系 ;在 
2-4 d 周 期 尺度 上 ,中 层 水 温 变 化 要 早 于 DO 浓度 变化 约 2~4h。 


关键 词 : SERÉ 冰 封 期 ; 水温 ; 溶解 氧 ; 小 波 变换 


溶解 氧 (DO ) 能 够 反映 水 体 污染 程度 ,是 水 生生 
物 生命 活动 不 可 缺少 的 环境 因子 ,同时 也 是 维系 水 
生生 物 存活 和 水 生态 环境 平衡 的 重要 因素 "'“”。 深 
解 氧 这 高 或 过 低 均 会 对 水 生生 物 尤 其 是 鱼 类 产生 
不 利 影响 。 当 水 体 中 溶解 氧 浓度 大 于 14.4 mg L`, 
导致 鱼 类 产生 泡 病 ,使 鱼 体 上 浮 或 游 动 失去 平衡 ; 
当 其 浓度 在 1~2 mg:L"' 时 ,大 部 分 鱼 类 的 生存 会 受 
到 威胁 ,但 当 其 浓度 变 高 后 , 鱼 类 可 恢复 健康 ; 若 其 
浓度 小 于 1 mg*L', 鱼 类 则 会 直接 死亡 , 底 泥 与 水 
体 之 间 的 物质 平衡 也 被 打破 ,有 害 物质 向 水 体 的 释 
放量 会 增加 。 

湖泊 DO 变化 过 程 受 温度 .水生 生物 呼吸 .水汽 
界面 的 气体 交换 、 光 合作 用 、 沉 积 物 及 底 栖 生 物耗 
氧 “ 等 多 种 因素 相互 影响 。 一 般 来 说 ,温度 的 变化 


0.4 mg"L'), 晶 温度 越 低 升温 后 溶解 度 下 降 越 明 显 ， 
尤其 对 于 一 些 寒带 和 温带 湖泊 而 言 ,水 温 上 升 引起 
的 DO 下 降 更 为 明显 ”另外 , 当 湖 泊 表 、 底 水 层 温 
度 相 差 较 大 时 ,表层 DO 很 难 扩散 至 湖泊 底部 ,再 加 
上 底部 有 机 物 浓度 较 高 , 耗 氧 快 ,使 得 底层 出 现 长 
时 间 的 缺 氧 和 厌 氧 状态 ”"。Jankowski 等 "基于 Zur- 
ich 湖 和 Greifensee W] 1956—2003 年 夏季 水 温和 DO 
剖面 测量 发 现 ,2003 年 湖 表层 与 底层 的 温差 比 
1956—2002 年 平均 值 高 出 2.2% 和 2.1% ,这 种 极端 
的 温度 垂直 梯度 极 大 地 阻碍 了 DO 的 扩散 和 上 下 水 
层 交 换 。 而 在 富 营养 化 浅水 湖泊 中 观测 发 现 , 水 温 
和 DO 分 层 的 现象 会 导致 内 源 营 养 盐 加 速 释放 ,使 
藻类 集中 于 水 体 表层 ,加 速水 华 暴发 "。 

寒 区 湖泊 具有 和 较 长 的 冰 封 期 , 冰 羡 阻隔 了 大 气 


会 直接 影响 水 体 中 DO 的 含量 ,特别 是 表层 DO 受 影 
响 特 别 明显 。Parmar 等 中 利用 回归 分 析 讨 论 了 在 印 
度 亚 称 纳 河 的 拦 河 坝 采样 点 监测 的 一 些 水 质 参 数 ， 
发 现 水 温和 深 解 氧 呈 负 相关 。 水 温 的 上 升 会 降低 
氧 在 水 中 的 洲 解 度 ,在 0~40 % 之 间 ,温度 每 上 升 
1 C, 氧 的 溶解 度 便 会 降低 1.6%~2.8%( 相 当 于 0.1~ 


收 稿 日 期 : 2020-07-22; 修订 日 期 : 2020-11-11 


与 湖水 的 气体 交换 ,强烈 的 垂直 混合 和 光合 作用 受 
到 抑制 ”。 水 的 通气 和 光合 作 用 产生 的 氧气 供应 
受到 限制 时 , 冰 下 水 体 中 氧气 就 会 逐渐 耗 光 , 造 成 
水 质变 化 并 产生 负面 的 生态 后 果 , 如 冬季 鱼 类 死 
亡 \ 底 枉 生物 丧失 以 及 厌 氧 过 程 的 激活 与 危险 化 合 
物 的 积累 。 在 冬季 水 体 DO 研 究 方面 ,Bai 等 "“ 
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通过 对 2011 年 1 月 13 日 至 4 月 25 日 的 Valkea- 
Kotinen 湖 冰 水 温和 DO 进行 Daubechies 小波 与 高 斯 
拟 合 发 现 ,在 4 月 5 日 前 ,水 柱 上 层 水 温 与 DO 浓度 
呈正 相关 关系 , 且 上 层 水 温 与 DO 的 相关 关系 更 接 
近 ; 从 融 冰 开始 到 冰 层 破裂 ,DO 浓度 随 水 温 升 高 而 
降低 。 李 明 等 中 用 Daubechies 小 波 对 寒 区 湖泊 各 深 
度 DO 进 行 多 尺度 分 析 , 通 过 分 析 DO 数 据 的 低频 和 
高 频 重 构 曲 线 来 研究 溶解 氧 的 变化 ,给 出 有 利于 数 
据 分 析 的 测量 时 间 间 隔 。 同 时 ,冬季 水 温 分 层 的 现 
象 会 逐渐 减弱 ,并 随 着 垂直 方向 的 水 温 梯度 变 小 ， 
DO0O 较 高 的 表层 水 体 下 沉 逐 渐 增 强 ,水 温 的 垂 向 差异 
逐渐 变 小 ”。 

乌 梁 素 海 地 处 北方 塞 旱 区 ,每 年 约 有 4~5 个 月 
的 冰 封 期 。 目 前 ,已 有 学 者 探究 了 该 水 域 冰 封 期 
DO 的 动态 分 布 与 温度 关系 。 孙 标 等 对 2009 年 乌 
梁 素 海 总 排 干 沟 过 饱和 溶解 氧 出 现 因素 进行 分 析 
发 现 ,DO 浓度 在 昼 间 随 着 光照 时 间 的 延伸 不 断 增 
加 , 且 DO 过 饱和 状态 易 出 现在 5 月 初 。 宋 更 "通过 
对 2015 一 2016 年 冬季 冰 下 水 体 DO 人 研究 发 现 ,DO 浓 
度 与 水 温 均 存 在 明显 的 分 层 现 象 和 昼夜 变化 模式 ， 
且 随 着 水 深 的 增加 ,DO 浓度 降低 。 但 温度 与 DO 之 
间 的 协调 变化 规律 与 机 制 还 尚 不 明确 ,因此 ,本 文 
就 乌 梁 素 海 冰 下 水 体 水 温 与 DO 浓度 变化 规律 进行 
研究 ,以 期 能 更 全 面 地 掌握 湖泊 溶解 氧 浓度 的 动态 
过 程 以 及 与 温度 的 关系 ,对 修复 富 营 养 化 水 体 有 重 
要 意义 Po21。 


1 材料 及 方法 


1.1 研究 区 概况 

EYER YS (40°36'~41°03'N, 108°43’~108°57’E) 
FEF RS ACA 7k C1) PFA Se T EE CK 
市 乌拉 特 前 旗 境 内 , 黄河 河套 平原 的 末端 , 西 临河 
套 灌区 , 东 靠 乌拉 山西 区 ,湖区 呈 南 北 狭长 .东西 束 
案 形 态 , 现 有 水 域 面积 325.31 km? ,其 中 芦苇 面积 ; 
202.2 km?, 明 水 区 面积 123.11 km?。 水 质 指 标 TN、 
TP、COD 等 一 度 达到 劣 V 类 ,湖泊 富 营 养 化 严重 。 
富 营 养 化 催生 了 以 龙 须 眼 子 菜 为 主 的 水 草 疯 长 ,水 
草 的 腐烂 会 使 营养 物 外 汇 , 湖 水 中 分 解 者 在 分 解 营 
养 物 质 时 会 消耗 大 量 溶解 氧 ,使 得 水 体 总 溶解 氧 下 
降 。 
1.2 数据 收集 与 处 理 

乌 梁 素 海 湖泊 L15 监测 点 位 于 乌 梁 素 海 湖区 中 
部 , 西 邻 茂盛 的 水 生生 物 , 回 东 靠近 湖岸 ,属于 湖水 
缓冲 区 , 且 水 深 较 浅 ,与 沉积 物 界 面 间 的 物质 交换 
频繁 .过程 复 杂 。 沉 积 物 的 耗 氧 约 占 整 个 水 体 耗 氧 
的 90% ,对 上 履 水 体 DO 含量 影响 较 大 。 当 水 体 长 
期 缺 氧 又 会 使 得 水 体 和 沉积 物 环 境 发 生 一 系列 变 
化 ,如 沉积 物 - 水 界面 宏观 底 栖 生物 消失 ,微生物 代 
谢 对 有 机 质 的 作用 从 好 氧 降解 变 为 厌 氧 降解 。 
此 选取 L15 为 本 次 试验 监测 点 ,于 2019 年 1 月 24 日 
0:06 分 至 2 月 24 日 23:56 对 监测 点 进行 长 达 32 d 的 
Kit DO 监测 (图 2), 该 阶段 属于 湖泊 冰 封 稳定 期 ， 
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。 城镇 
-一 排 干 与 水 系 


图 1 研究 区 概况 示意 图 


Fig. 1 Overview of the research area 
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一 排 干 
= 明 水 区 
mm 芦苇 区 


图 2 采样 点 位 置 示意 图 


Fig. 2 Location diagram of sampling points 


湖面 冰 层 相对 稳定 。 同 时 ,在 2 月 13 一 17 日 有 两 场 
降雪 ,第 一 场 降雪 的 最 大 雪 厚 约 为 S cm, 第 二 次 降 
雪 的 最 大 雪 厚 约 为 3 cm( 表 1)。 

依据 监测 点 实际 水 深 ,划分 为 表 、 中 、 底 等 3 
层 。 其 中 ,将 冰 下 20 cm 的 表层 水 体 划 分 为 表层 ; 湖 
水 - 底 泥 交 界 处 以 上 约 10 cm 定 为 底层 ;中 层 则 根据 
监测 点 实际 水 深 的 二 分 之 一 。 本 次 测量 分 别 在 表 、 
中 、 底 3 个 深度 悬挂 TDO 探头 (深度 的 零点 参考 水 
平 是 冰 层 下 履 水 体 ) ,监测 时 间 间 隔 为 10 min。 采 用 
Morlet 小 波 分 析 法 分 析 水 温 与 DO 浓度 变化 特征 的 
周期 性 ,根据 小 波 系数 的 实 部 .小波 方差 小波 相干 
谱 等 信息 揭示 乌 梁 素 海 冰 封 期 冰 下 表层 水 体 (0 
em) 中 层 水 体 (66 em) .底层 水 体 (130 cm) 的 水 温 
和 溶解 氧 (DO ) 浓 度 的 特征 规律 。 


2 小 波 变 换 原 理 


小 波 分 析 广 泛 应 用 于 多 时 间 尺 度 变 化 特征 的 
研究 中 ,其 目的 是 得 到 发 生 在 不 同时 间 尺 度 下 的 局 
部 和 瞬 态 现象 的 完整 表示 站。 小 波 变 换 (Contin- 
ue Wavelet Transform ) 包 括 连 续 小 波 变 换 和 离散 小 
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波 变换 。 本 文 所 分 析 的 数据 是 水 体温 度 与 溶解 氧 
浓度 的 时 间 序 列 数据 ,实际 工作 中 ,该 时 间 序 列 的 
数据 常常 是 离散 的 ,因此 小 波 变 换 (Continue Wave- 
let Transform) KY A HOB SUE LA: 
W (a. b) -|aE [7 f oeh ^r (1) 

式 中 :a 为 尺度 因子 ;5 为 时 间 因 子 ; W,(a, b) 为 变换 
系数 。 

小 波 方差 定义 ,将 时 间 域 上 的 关于 a 的 所 有 小 
波 系数 的 平方 进行 积分 , 即 为 小 波 方差 , 它 能 反映 
言 号 波动 能 量 随 尺 度 a 的 分 布下。 小 波 方 差 图 可 用 
来 确定 信号 中 不 同 尺 度 扰动 的 相对 强度 和 存在 的 
主要 时 间 尺 度 , 即 主 周期 。 

V,(a)= |W (a, bf db (2) 

小 波 相干 谱 可 以 用 来 度量 时 频 空 间 中 两 个 时 
间 序 列 局 部 相关 的 密切 程度 ” ,其 原理 是 通过 对 两 
个 时 间 序 列 采 用 交叉 小 波 变化 (Cross Wavelet Sec- 
tumy) 的 方法 进行 连续 小 波 变 换 (CWT) ,以 此 来 揭示 
这 两 个 时 间 序 列 共 同 的 高 能 量 区 域 。 其 中 ,时 频 空 
间 中 的 两 个 时 间 序 列 局 部 相关 程序 需要 小 波 相 干 
谱 (Wavelet Coherence ) {jo X 时间 序 列 和 YY 时 间 
序列 的 小 波 相干 谱 函 数 为 : 


| S( wp (2) | 
s" W; (s) j x SG W, (s) b 
式 中 :3$ 是 平滑 器 ;* 为 伸缩 尺度 ; Wz (s). W, (s) 分 别 
为 XY 的 小 波 变 换 ; WW (s) 为 交叉 小 波谱 。 


3 ARAM 


3.1 温度 溶解 氧 的 分 布 特征 

依据 所 监测 的 原始 数据 绘制 温度 、 溶 解 氧 的 
垂 向 分 布 趋势 图 (图 3a、 图 3b)。 由 图 3 可 以 看 出 ， 
温度 有 明显 的 层 化 结构 , 即 表 -中 -底层 , 冰 下 表层 
水 体 平均 水 温 达 -0.12 ;中 层 平均 水 温 达 
5.79 % ;底层 平均 水 温 达 6.57 °C ;水 体 整体 温度 浮 
动 范围 在 -0.13~7.66 包 。 表 层 最 高 水 温 出 现在 2 月 


Ri(s)= 


表 1 冬季 野外 观测 记录 表 


Tab. 1 Winter field observation record table 


年 份 冰 厚 /m 最 大 雪 厚 /cm 降雪 覆盖 期 多 云 期 
2 月 13 一 15 日 2 月 4 一 10 日 
ae ae us ? 2 月 17 一 20 日 2 月 13 一 24 日 


202106.00060v1 


chinaXiv 


632 T 


旱 


Bl 


研 


ChinaXiv 合 作 期 刊 


38 卷 


究 


时 间 
01-240:06 102-030:06 02-130:06 02-230:06 ÜC 


m 7.660 
6.753 
5.845 
4.938 
4.030 
3.123 
2.215 
1.308 
0.4000 


深度 /cm 


深度 /cm 


时 间 
01-24 0:06 02-03 0:06 02-13 0:06 02-23 0:06 DO/(mgL?) 


(b) 12.95 
20 11.33 
40 9.713 
8.094 

60 
6.475 
80 4.856 
100 3.238 
1.619 

120 
0.000 


图 3 温度 (a) 溶解 氧 (b ) 剖 面 结构 


Fig. 3 Temperature (a), dissolved oxygen (b) section structure diagram 


24 H 13:06(0.34 © ) ,最 低 水 温 出 现在 1 月 24 日 
00: 26 ( -0.13 C ); 中 层 最 高 水 温 出 现在 2 月 7 日 
11:36(7.37 C) ,最 低 水 温 出 现在 1 月 25 日 6:46 
(3.58 °C ); 底 层 最 高 温度 出 现在 2 月 7 日 15:16 与 
15:36(7.66 °C) ,最 低 水 温 出 现在 1 月 24 日 00:56 
(5.47% )。 且 三 层 水 体温 度 存在 着 底层 温度 偏 高 、 
表层 温度 偏 低 的 分 层 趋势 。 

溶解 氧 剖 面 (图 3b) 与 水 温 剖 面 的 分 层 不 一 致 ， 
DO 浓度 在 不 同 水 深 下 的 分 布 情况 不 同 。 表 层 DO 
浓度 浮动 在 9.12~12.11 mg*L', 平 均 浓度 达 10.91 
mg:L', "HZ DO YKE JJ TE 0.90-13.44 mg*L', 平 
均 浓度 达 7.14 mg: L'o JEJ DO 浓度 浮动 在 0.04~ 
4.72 mg: L'' 平均 浓 度 达 1.82 mg.L'。 由 图 3 可 知 ， 
上 混合 层 、 中 间 层 、 下 混合 层 也 存在 明显 分 层 现 
象 。 表 层 DO 浓 度 近 1 个 月 时 间 内 基本 保持 不 变 ;而 
中 层 DO 浓 度 从 2 月 4 日 8:26 开 始 逐 渐 下 降 , 由 
10.18 mg:L' FEES 3.82 mg L'', FA E AR 0.32 
mg:L'«d', HITS SIE TE us e C UTERIS , 
冰 下 水 体 对 光 的 透射 率 下 降 到 1% 左 右 ,湖泊 水 
体 的 光合 作用 减弱 , 产 氧 量 低 ,有 机 物 与 生物 的 呼吸 
作用 占 主导 地 位 ,DO 被 大 量 消 耗 ,中 层 氧 含量 发 生 反 
转 ;底层 DO 浓度 由 1.93 mg FIAESI 1.11 mgsL^, F 
降 速率 为 0.05 mg*L!'.d'。 表 层 溶 解 氧 浓度 浮动 幅 
度 较 平 绥 , 随 着 冰 封 期 的 时 间 逐 渐 增 长 ,中 、 底 层 深 
解 氧 下 降幅 度 相 对 明显 。 且 水 体 中 的 DO 浓度 也 出 
现 了 分 层 现象 。 
32 水 温 与 DO 浓度 小 波 变换 分 析 
32. 水 温 的 小 波 变换 分 析 ”为 进一步 揭示 乌 梁 素 
海 冰 封 期 水 体 分 层 的 温度 在 时 间 尺 度 上 变化 的 多 
样 性 。 采 用 Morlet 小 波 进行 32 d 尺 度 的 一 维 连续 小 


波 变换 获得 小 波 系数 实 部 等 值 线 与 小 波 方 差 图 进 
行 多 时 间 尺 度 变化 分 析 。 由 表层 水 温 Morlet 小 波 实 
部 等 值 线 ( 图 4) 可 知 ,表层 水 体温 度 演变 过 程 存在 
多 时 间 尺 度 特征 。 其 中 相对 比较 明显 的 特征 时 间 
尺度 为 9~17 d、19~29 d, 两 类 尺度 的 周期 变化 规 
律 。19~29 d 特征 时 间 尺 度 在 2 月 22 日 之 前 存在 
高 - 低 - 高 - 低 的 3 次 震荡 ,随后 呈 增 强 态 势 ;9~17 d 
尺度 以 16 d 为 尺度 中 心 ,ARERE F 19~29 d 
度 周 期 变化 中 。 另 外 ,通过 表层 水 体温 度 的 小 波 方 
差 大 小 的 波动 ,揭示 了 水 温 的 周期 震荡 强度 。 小 波 
方差 存在 2 个 较为 明显 的 峰值 ,依次 对 应 着 16 d、25 
d 时 间 尺 度 。 其 中 最 大 峰值 对 应 25 d 的 时 间 尺 度 ， 
说 明 25 4 左右 的 周期 震荡 最 强 ,为 表层 水 体温 度 变 
化 的 第 一 主 周 期 ;16 d 时 间 尺 度 对 应 着 第 二 峰值 ,为 
水 体温 度 变 化 的 第 二 主 周期 ,2 个 周期 的 波动 控制 
着 表层 水 体温 度 在 整个 时 间 域 内 的 变化 特征 ,从 侧 
面 验证 了 实 部 等 值 线 图 中 振荡 周期 。 

中 层 的 水 体温 度 变 化 存在 多 尺度 的 周期 变化 
规律 ,其 中 演变 过 程 较 明显 的 特征 时 间 尺 度 周期 是 
10-17 d, 在 该 特征 时 间 尺 度 上 ,中 心 时 间 尺 度 为 14 
d, 水 温 变 化 明显 ,频率 快 ,水 温 交 替 变 化 ,具有 全 局 
性 ;而 在 25~32 d 特 征 时 间 尺 度 上 水 温 呈 现 的 是 低 - 
高 - 低 的 变化 ,能 量 波及 大 ,但 未 形成 闭合 圆圈 ,不 
能 作为 特征 周期 进行 分 析 。 小 波 方差 大 小 的 波动 
揭示 了 中 层 水 温 的 周期 震荡 强度 。 在 10~17 d 时 间 
尺度 下 水 温 浮动 较 大 ,中 心 时 间 尺 度 14 d 是 表层 温 
度 变化 的 第 一 主 周期 。 

底层 水 温 小 波 实 部 图 等 值 线 显示 底层 水 温存 
在 两 类 尺度 的 周期 变化 规律 。1~6 4 尺度 ,以 2.5d 
为 震荡 中 心 ,水 温 变化 强度 较 弱 ,震荡 能 量 平稳 ;9~ 
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图 4 水 体温 度 的 小 波 系 数 实 部 等 值 线 图 及 小 波 方 差 


Fig.4 Realpart contour map of wavelet coefficient and wavelet square of water temperature 


19 d 尺 度 ,作为 32 d 内 主要 周期 变化 ,震荡 中 心 为 
14 d, 出 现 了 明显 的 水 温 升 降 交 蔡 周期 ,日 在 该 周期 
下 存在 一 个 高 温 期 。 高 能 量 中 心 集中 在 2 月 13 一 16 
日 之 间 , 时 间 跨 度 约 3 d ,边缘 闭合 程度 强 ,说明 该 尺 
度 下 水 体温 度 有 规律 的 增 大 ,通过 小 波 方差 波动 
图 ,进一步 揭示 温度 变化 的 周期 震荡 强度 。 小 波 方 
差 中 存在 2.5 d、14 d 两 个 峰值 ,14 d 时 间 尺 度 为 最 大 
峰值 ,是 水 温 变化 第 一 主 周 期 ,第 二 主 周期 则 为 2.5 
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时 间 尺 度 /d 


d, 第 二 主 周期 能 量 波动 小 ,第 一 主 周期 能 量 波动 
大 ;从 侧面 验证 了 实 部 等 值 线 图 的 震荡 周期 。 
根据 小 波 方 差 检 验 的 结果 ,绘制 水 体 各 层 水 温 
的 第 一 主 周期 小 波 系 数 图 (图 5)。 从 主 周期 趋势 图 
中 我 们 发 现在 特征 时 间 尺 度 下 ,表层 水 体 水 温存 在 
升 高 -降低 变化 特征 。 在 25 d 特征 时 间 尺 度 上 ,水 
温 变化 的 平均 周期 为 15 d 左 右 , 大 约 经 历 了 两 个 
升 - 降 转换 期 ;在 小 波 变 换 中 ,不 同时 间 尺 度 小 波 系 


图 5 水 体温 度 的 特征 时 间 尺 度 小 波 实 部 过 程 线 


Fig. 5 Characteristic time scale wavelet real part process line of water body temperature 
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数 ( 实 部 ) 过 程 线 的 过 零点 对 应 突变 点 ,其 中 低频 尺 
度 更 能 反映 序列 实际 变化 。 中 层 水 体 水 温 在 14 d 
特征 时 间 尺 度 下 ,3 个 升 - 降 周期 变化 ,水 温 变 化 的 
周期 约 为 9 d; 底 层 水 体 水 温 在 第 一 主 周期 14d 尺 度 
上 发 生 3 次 突变 :1 月 28 日 ( 升 - 降 )、2 月 7 日 ( 升 - 
降 )、2 月 15 日 ( 升 - 降 ), 其 变化 周期 约 为 8.5 do 
3.2.2 DO 浓度 的 小 波 变 换 分 析 DO 浓度 的 小 波 系 
数 实 部 等 值 线 图 反映 了 序列 在 不 同时 间 尺 度 的 周 
期 变化 、 时 域 分 布 ,以 及 时 间 序 列 的 变化 趋势 (图 
6) ,颜色 越 蓝 表示 DO 浓度 越 低 。 颜 色 越 浅 则 该 时 
期 DO 浓度 高 , 零 值 对 应 突变 点 。 图 6 可 以 发 现 : 表 
层 水 体 DO 浓 度 存在 1 类 周期 变化 ,19~26 d 时 间 尺 
度 , 以 23 d 为 震荡 中 心 ,振荡 能 量 偏 弱 。 

通过 中 层 DO 浓度 等 值 线 图 可 以 看 到 ,在 6~11 
d 时 间 尺 度 上 ,DO 浓度 的 变化 程度 不 高 , 仅 小 幅 变 
动 ; 而 在 15~23 d 时 间 尺 度 上 ,DO 浓度 从 1 月 24 日 


至 2 月 24 日 的 32 d 内 ,有 3 个 高 浓度 时 期 和 3 个 低 
浓度 时 期 ,1 月 24 日 前 后 进入 第 一 个 低 浓 度 期 ,2 月 
18 日 前 后 进入 最 后 一 个 低 浓度 期 。 该 周期 变化 下 ， 
中 层 水 体 D0 浓 度 变化 呈现 降低 -上 升 -降低 -上 升 - 
降低 -上 升 的 变化 过 程 。 在 30 d 时 间 尺 度 上 ,DO 浓 
度 出 现 高 低 交 蔡 的 3 次 震荡 ,该 周期 变化 主要 活跃 
于 1 月 24 日 至 2 月 22 日。 为 得 到 中 层 水 体 DO 浓 度 
值 在 时 间 序 列 的 主 周期 情况 , 画 出 一 维 连 续 小 波 变 
换 系数 相应 的 小 波 方 差 图 。19 d 尺 度 上 有 一 个 比较 
强 的 波动 能 量 ,在 7.5 d 尺 度 左 右 有 一 个 较 强 的 波动 
能 量 。 而 时 间 尺 度 小 于 5 d 的 尺度 上 波动 能 量 相对 
而 言 很 小 。 这 表明 , 近 32 d 里 ,中 层 水 体 DO 值 演变 
在 该 尺度 上 ,存在 一 个 19 4 变换 的 第 一 主 周期 和 
7.5 d 变 化 的 第 二 主 周期 。 

对 底层 水 体 DO 浓度 进行 Morlet 小 波 分 析 , 得 出 
乌 梁 素 海底 层 水 体 32 d DO 浓度 值 的 变化 情况 。 在 
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图 6 水 体 DO 浓 度 的 小 波 系数 实 部 等 值 线 图 及 小 波 方差 


Fig.6 Real-part contour map and wavelet square of DO concentration in water body 
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19~29 d 时 间 尺 度 上 ,有 3 个 明显 的 高 浓度 时 期 。 但 
在 1 月 29 日 与 2 月 14 日 水 体 DO 开 始 进 入 低 浓 度 时 
期 。 尺 度 中 心 为 24 d, 在 该 时 间 尺 度 下 DO 浓度 变 
化 呈现 交替 周期 规律 ,处理 该 尺度 上 的 小 波 方差 
图 ,底层 水 体 DO 浓 度 值 在 24d 尺度 上 有 强烈 的 能 
量 波动 。 而 在 10 d 尺 度 上 ,能 量 波动 弱 , 变 化 不 明 
显 。 根 据 底层 DO 的 小 波 方差 图 ,可 以 看 出 底层 水 
体 DO 浓 度 值 存在 2 个 波峰 。24 d 时 间 尺 度 为 最 大 
峰值 ,是 DO 浓度 值 变化 的 第 一 主 周期 ;而 10 d 则 是 
DO 浓度 值 变化 的 第 二 主 周期 。 在 整个 时 间 域 内 ， 
底层 水 体 DO 浓 度 值 变化 主要 受到 上 述 两 种 周期 
制约 。 

根据 小 波 方差 检验 的 结果 ,绘制 水 体 各 层 DO 
浓度 变化 的 第 一 主 周期 小 波 系数 图 及 升 高 -降低 变 
化 特征 (图 7)。 表 层 DO 浓度 在 23 d 特 征 时 间 尺 度 
上 ,变化 的 平均 周期 为 15 d 左 右 ,大 约 2 个 周期 的 
升 - 降 变化 ;中 层 DO 浓度 在 8d 特 征 时 间 尺 度 下 ,经 
历 约 6 个 升 - 降 变 化 ,平均 周期 为 5 d; 从 底层 水 体 
DO 浓度 主 周期 趋势 图 中 我 们 发 现在 特征 时 间 尺 度 
24 d 的 变化 下 ,DO 存在 2 次 突变 ,分 别 是 :1 月 24 日 
至 2 月 1 日 发 生 浓度 的 升 高 -降低 突变 ,2 月 9 一 16 日 
发 生 第 二 次 突变 。 同 时 ,底层 DO 浓度 变化 的 平均 
周期 为 7d 左 右 。 
3.3 水 温 与 DO 浓度 的 小 波 相干 谱 分 析 

选取 了 2019 年 1 月 24 日 至 2 月 24 日 乌 梁 素 海 
的 L15 点 冰 下 表层 中层、 低层 水 体温 度 和 DO 浓度 
值 数据 ,利用 小 波 相干 谱 绘 制 温度 -DO 频谱 图 。 小 
波 相干 谱 揭 示 两 组 信号 序列 的 时 频 相 关 性 ,特别 适 
用 于 2 种 过 程 具有 很 强 的 相互 作用 时 的 时 间 序 列 时 
频 分 析 。 图 8 是 1 月 24 日 至 2 月 24 日 的 水 温 序列 和 
DO 浓度 序列 显示 的 在 2~8 d 的 周期 频段 内 的 小 波 
相干 谱 。 

图 8a 中 显示 :表层 水 体温 度 与 DO 浓度 序列 位 
相 在 2 月 15( 图 中 第 22 d) 一 24 日 期 间 , 在 3~7 d 的 周 
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图 7 DO 浓度 的 特征 时 间 
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期 区 间 表 现 出 通过 95% 的 置信 和 度 检验 的 显著 共振 
关系 (图 中 黑 粗 线圈 部 分 ) ,并 且 在 周期 区 间 内 箭头 
指向 基本 一 致 (向 左 ) ,表明 了 序列 的 反 相位 变化 ， 
相关 系数 达到 0.9; 图 8b 显示 出 ,中层 水 温 与 DO TK 
度 两 序列 在 2~4 4 尺度 的 显著 相关 性 的 特征 : 即 在 2 
月 9 日 (图 中 第 15 d) 开 始 呈 现 出 很 强 的 相关 关系 ， 
相关 系数 达到 0.9, 通 过 95% 的 置信 和 度 检 验 , 随 着 冰 
封 期 的 推进 ,水温 位 相 系统 开始 超前 DO 浓度 变 
化 。 经 过 计算 水 温 与 DO 浓度 在 2 月 11( 图 中 第 18 
d) 一 18 日 (图 中 第 25 d) 时 段 及 2~4d 的 日 变化 尺度 
上 ,水 温 序 列 变化 超前 DO 浓度 变化 约 2~4h; 图 8c 
显示 ,1~3 d 周 期 区 间 , 底 层 两 序列 在 2 月 20( 图 中 第 
28 d) 一 24 日 上 通过 95% 置 信和 度 的 检验 ,相关 系数 约 
为 0.8 ,箭头 的 方向 基本 上 在 0~45$。 ,表明 底层 水 温 
与 DO 浓度 在 2 月 20—24 日 上 具有 很 强 的 负 相关 关 
系 。 整 体 来 看 , 表 、 中 、 底 三 层 的 水 体温 度 与 DO 浓 
度 变 化 呈现 不 同 规律 。 


4 讨论 


对 于 中 国 北方 寒 旱 区 湖泊 来 说 ,有 无 冰 封 期 的 
DO 差异 性 主要 体现 在 : 湖 冰 下 的 热 分 层 颠倒 ,水 温 
较 低 ,空气 -水 -气体 交换 、 风 的 混合 作用 和 大 气 输 
入 基本 停止 。 通 过 对 湖泊 冰 下 水 体 的 DO 含量 与 温 
度 监 测 发 现 :底层 水 温 最 高 ,DO 含量 最 低 ; 表 层 水 
温 最 低 、.DO 含 量 最 高 。 通 常 来 讲 , 表 层 温 度 的 下 降 
是 由 于 冰 的 热量 损失 ,而 底层 温度 的 增加 则 是 来 自 
沉积 物 的 热 通 量变 化 "。 由 此 可 知 , 当 水 体 的 初始 
温度 越 低 , 冰 下 可 储存 的 DO 含量 就 越 高 。DO 垂 直 
方向 的 分 布 规律 则 是 表层 含量 高 .底层 含量 低 ,在 
分 层 湖 泊 中 底层 DO 含量 低 的 现象 是 很 普遍 的 ,例如 
美国 北部 的 Erie 湖 夏季 湖底 缺 氧 可 追溯 千年 ”。 冰 
封 期 内 半数 的 DO 消耗 发 生 在 冰期 的 第 一 个 月 内 ， 
表 - 中 - 底 的 DO 含量 呈现 出 由 表 及 底 逐 渐 下 降 的 趋 


尺度 小 波 实 部 过 程 线 


Fig.7 Characteristic time scale wavelet real part process line of DO concentration 
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周期 /a 


5 10 15 20 25 30 
时 间 /d 
(b) 温度 -DO 


5 10 15 20 25 30 
时 间 /d 
(c) 温度 -DO 


周期 /a 


5 10 15 20 25 30 
时 间 /d 


注 :箭头 表示 两 者 之 间 的 相位 关系 ,一 代表 两 个 时 间 序 列 之 间 的 反 相 位 关系 , 即 为 负 相关 ;一 代笔 两 个 时 间 序 列 
之 间 是 同 相 位 关系 , 即 为 正 相关 ; T 表示 后 者 变化 超前 者 变换 90°; | 表示 后 者 变化 落后 前 者 90°。 
粗 实 线 区 域 表示 显 著 性 a=0.05 条件 下 的 红 噪 声 标准 谱 的 检验 ,小 波 影 响 椎 ( 细 弧 线 ) 以 内 区 域 为 有 效 谱 值 。 
图 8 表层 (a)、 中 层 (b)、 底 层 (c) 的 水 温 与 DO 浓度 变化 的 小 波 相 干 谱 


Fig. 8 Wavelet coherence spectrum of water temperature and DO concentration changes in 


surface layer (a), middle layer (b) and bottom layer (c) 


PRAT AE TE RD AE th], DOR TE AYE To] EEA IL 
受到 生物 化 学 的 控制 ,更 重要 的 是 水 温 的 影响 ”， 


势 , 表 层 水 体 受 光合 作用 的 影响 ,产生 氧气 和 有 机 
物 ,使 得 该 层 DO 含量 升 高 ,但 部 分 有 机 物 会 沉降 至 


湖泊 底部 ,在 其 下 降 过 程 中 不 断 的 被 分 解 ,消耗 掉 
底层 水 体 D0 ,使 底层 DO 含量 降低 。 与 此 同时 ,无 
雪 的 冰 是 冬季 中 进行 光合 作用 的 先决 条 件 ,但 是 只 
要 有 效 光 合 辐射 可 用 ,DO 就 会 产生 ™。 因 此 ,我 们 
认为 光合 作用 是 受 太阳 辐射 的 控制 ,也 受到 冰 层 上 


体现 在 表层 温度 低 (-0.12 °C) EB dl ES 
(6.57 C), 表 、 底 的 温差 阻碍 了 水 体 对 流 ™ 进 而 波及 
到 中 层 DO 浓 度 。 此 外 ,湖泊 又 受到 气候 .地理 、 富 
营养 化 有 机 污染 物 .水 文 及 负荷 的 输入 过 程 等 环 
境 压 力 的 作用 ,使 分 层 后 的 DO 在 表 、 中 和 底层 的 浓 
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度 形态 与 含量 发 生变 化 。 大 气 复 氧 .光合 作用 等 会 
导致 表层 DO 浓度 较 大 ,但 随 着 冰 封 期 冰 层 的 增 
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(2) 乌 梁 素 海 水 体温 度 与 DO 浓度 变化 周期 变 
化 呈现 高 低 交 替 的 位 相 特 征 ,小 尺度 高 低位 相 变 化 


长 以 及 冰雪 覆盖 ,阻隔 大 气 与 水 体 的 气体 交换 , 同 
时 限制 表层 水 体 水 生生 物 的 光合 作用 ,造成 在 一 
定时 期 内 DO 含量 下 降 。 而 底层 的 缺 氧 现象 则 是 与 
沉积 物 、 初 级 生产 力 、 外 部 输入 和 克 (P) 的 释放 相 
Ko Zhao 等 "在 乌 梁 素 海中 观察 到 的 沉积 物 ,他 们 
推测 ,沉积 物 中 释放 出 4 mg*L' 的 P, 表 明 沉 积 物 内 
积累 的 磷 含 量 可 能 高 于 冬季 水 柱 中 的 平均 磅 含 
量 。 同 时 ,湖泊 分 层 后 底层 通常 表现 为 耗 氧 状态 。 

正 是 由 于 冰 下 水 体 水 温存 在 分 层 ,导致 D0 含 
量 也 出 现 分 层 , 水 温 增加 会 使 DO 浓度 降低 ,水 温 的 
降低 亦 会 使 DO 浓度 增加 。 这 一 规律 说 明了 水 温 对 
于 水 体 中 的 DO 浓度 的 影响 ,是 一 种 负 相 关 关系 的 
绝对 影响 ,通过 小 波 相 干 谱 的 结果 就 可 以 看 到 , 表 、 
底 的 水 体温 度 和 DO 浓度 的 相位 关系 是 负 相 位 关 
系 , 相 关系 数 达 0.9, 表 明 两 者 之 间 存 在 高 强度 的 负 
相关 关系 。 并 且 在 水 温 的 控制 下 物理 和 生化 过 程 
都 会 对 DO 浓度 造成 影响 。 水 温 低 时 ,会 抑制 水 中 
有 机 物 与 生物 的 呼吸 作用 。 使 DO 含量 保持 一 定 程 
度 的 稳定 ; 而 沉积 物 -水 界面 的 温度 较 高 ,接近 乌 深 
素 海 地 区 多 年 的 平均 气温 ,温度 的 适宜 性 大 大 文 持 
了 沉积 物 及 底 枉 生 物 的 耗 氧 ,使 DO 浓度 下 降 ,底层 
逐渐 变 成 缺 氧 层 。 减 少 的 DO 会 通过 扩散 作用 去 影 
响 整 个 水 柱 :使 中 层 水 温 与 DO 浓度 由 原来 的 负 相 
关 关 系 逐 渐 转 变 为 水 温 超前 DO 浓度 变化 2~4 ho 
且 随 着 冰 封 时 间 的 持续 ,DO 浓度 会 越 降 越 低 , 冰 下 
水 体 中 缺 氧 层 离 湖底 的 距离 也 越 来 越 远 。 

文中 仅 基于 Morlet 小 波 分析 了 冰 封 期 乌 深 素 海 
分 层 水 体温 度 与 DO 浓度 的 演变 规律 ,但 关于 冰 下 
水 体温 度 与 DO 浓度 之 间 的 深层 耦合 机 理 还 有 待 我 
们 利用 更 长 时 段 的 数据 来 进一步 分 析 和 研究 。 


5 结论 

(1) 乌 梁 素 海 冰 下 水 体 层 按 水 温 DO 浓度 变 
化 趋势 可 分 为 表 .中 、 低 三 层 , 表 层 水 体 平均 水 温 达 
-0.12 % ;中 层 平均 水 温 达 5.79 % ;底层 平均 水 温 达 
6.57 CC。 表层 水 体 DO 浓 度 平均 值 为 10.91 mg-L"'; 
中 层 DO 浓度 平均 值 为 7.14 mg: L'; JEJE DO 浓度 平 
均值 为 1.82 mg:L'。 自 2 月 4 日 ,DO 的 消耗 增 大 。 
水 体 DO 浓 度 含量 亏损 未 得 到 补充 ,整个 水 柱 层 的 
DO 含量 发 生 改 变 。 


舱 套 于 大 尺度 位 相 变 化 中 。 表 层 水 温 变化 的 主 周 
期 是 25 d.16 d, 在 其 第 一 主 周 期 下 ,大 约 经 历 了 2 个 
升 - 降 转 换 期 ,周期 约 为 15 d; DO 浓度 变化 的 主 周期 
是 23 d, 同 样 经 历 了 2 个 周期 的 升 - 降 变化 ,平均 周 
期 为 15 d 左 右 ;中 层 水 温 变 化 的 主 周期 是 14 d, DO 
浓度 变化 的 主 周 期 是 19 d、7.5 d, 在 第 一 主 周期 下 ， 
水 温 3 个 升 - 降 周期 变化 ,平均 周期 约 为 9d,DO 深 
度 经 历 约 6 个 升 - 降 变化 ,平均 周期 为 5 d; 底 层 水 温 
变化 的 主 周期 是 14 d、2.5 d, DO 浓度 变化 是 24 d, 
10 d。 底 层 水 体 水 温 在 第 一 主 周期 下 ,水 温 发 生 3 
次 突变 ,变化 周期 约 为 8.5 d, DO 浓度 发 生 2 次 突变 ， 
平均 周期 7d 左 右 。 除 了 底层 变化 周期 外 ,其 他 TT- 
DO 变化 周期 大 臻 相近。 水 温 高 能 量 区 的 分 布 几乎 
也 是 DO 浓度 高 能 量 区 的 集合 ,水 温 的 变化 对 DO 的 
变化 有 着 积极 的 响应 。 

(3) 温度 -DO 小 波 相干 谱 分 析 显 示 ,整个 水 柱 
三 层 的 水 温 .DO 浓度 周期 变化 不 相同 。 同 时 , 冰 封 
期 乌 梁 素 海 各 水 层 的 温度 与 DO 浓度 之 间 存 在 复杂 
的 关系 :表层 水 体温 度 与 DO 浓度 序列 位 相 在 3~7 d 
的 周 反 相位 变化 ,相关 系数 达到 0.9; 中 层 水 温 与 DO 
浓度 两 序列 在 2~4 d 尺度 在 2 月 9 日 开始 呈现 出 很 
强 的 负 相 关 关 系 , 随 着 冰 封 期 的 推进 ,水 温 位 相 系 
统 开 始 超 前 DO 浓度 变化 约 2~4 h; 1-3 d 周 期 区 间 ， 
底层 两 序列 通过 95% 置 信 度 的 检验 ,相关 系数 约 为 
0.8 ,底层 水 温 与 DO 浓度 存在 负 相 关 关 系 。 
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Change law of water temperature and dissolved oxygen concentration of 
Wuliangsu Sea in icebound period 


ZHAI Jialun, SHI Xiaohong, LIU Yu, ZHAO Shengnan, BAO Wenzhi LI Guohua 
(College of Water Conservancy and Civil Engineering, Inner Mongolia Agricultural University, 
Hohhot 010018, Inner Mongolia, China) 


Abstract: The wavelet analysis method was used to analyze the trend of temperature and oxygen concentration 
of a water body under ice in Ulansuhai in 2019, and the time series was one month: from January 24 to February 
24, 2019. The results showed that the water temperature and oxygen concentration of the lake were stratified. The 
oxygen concentration in the surface layer did not change. However, the oxygen concentration decreased in the 
middle (0.32 mg: L~- d") and bottom (0.05 mg-L d) layer. The main period of static water temperature in the 
surface layer was 25 d and 16 d. Under the first principal period, the period was about 15 d. The main period of 
static oxygen concentration was 23 d, and the average period was 15 d. The main period of static water 
temperature in middle layer was 14 d, and the main period of static oxygen concentration was 19 d and 7.5 d. 
Under the first principal period, the average period of static water temperature was about 9 d, and the average 
period of static oxygen concentration was 5 d. The main period of water temperature in the bottom layer was 14 d 
and 2.5 d, and the oxygen concentration was 24 d and 10 d. The period of change in water temperature was about 
8.5 d, and the average period of oxygen concentration was about 7 d. The wavelet coherence spectrum results 
showed that there is a complex relationship between lake water temperature and the oxygen concentration, 
suggesting the water temperature shows an inverse phase relationship with the oxygen concentration in the 
surface and bottom layer of the lake. Additionally, on the 2-4 d cycle scale, the change of water temperature in 
the middle layer was about 2-4 h ahead of the change in dissolved oxygen concentration. 


Keywords: Wuliangsu Sea; icebound period; water temperature; dissolved oxygen; wavelet transform 


